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ABSTRACT 

Surface wave action has been proven to be the biggest deterrent 
to a ship's progress through the water. Until the present time fore- 
casts of surface waves have been available for only a few days. In 
order to determine an optimum ship route across the ocean these wave 
forecasts must be extended to 15-20 days. Using three computer programs 
from the Fleet Numerical Weather Facility and the U.S. Weather Bureau 
5-day and 30-day forecasts of the mean sea level pressures as input 
data, surface wave forecasts may be prepared for periods up to 30 days. 
When compared with the actual waves over the specified period of time, 
the forecast wave patterns show more skill than forecast waves derived 


strictly from climatology. 
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1. INTRODUCTION 

Optimum track ship routing (OTSR) has been in existance in the 
U.S. Navy for almost ten years under this specific name. An optimum 
track 1S one chosen from port to port across the ocean that will allow 
a ship to traverse the route in the least amount of time and also en- 
counter less unfavorable weather. James [5 | has shown that a definite 
increase in operational efficiency is made possible by the routing of 
transoceanic vessels along tracks based primarily on long-range ocean 
wave forecasts rather than along standard seasonal tracks. Such predic-~ 
ted tracks are believed to be superior to the standard seasonal routes 
in terms of travel time, passenger comfort, damage to cargo and ship, 
and fuel consumption. Absolute values for the reduction of the above 
items are difficult to present, but certainly the great weight of evi-~ 
dence presented in his evaluation supports the contention that efficiency 
of ship operation is improved by oceanographic ship routing. The Mili- 
tary Sea Transportation service tentatively estimates that each hour 
saved in ship operations is a savings of $200.00 to them. 

Some of the international steamer tracks in use today are modifi- 
cations of routes laid out by Lt. Matthew Fontaine Maury, USN, as far 
back as the 1850's. Maury's routes were based on ships' logs, clima-~ 
tology and statistics. Today, james [4} has found that the most impor- 
tant parameter retarding a ship's progress is Surface wave action. 

Wave action actually affects a ship's progress in two ways. First 
we by retarding the ship's speed through the disturbed water and 
secondly by making the Captain order a reduction in speed due to the 
violent effects of roll, pitch, and other associated motion of the ship 


in creating eddy resistance. 


Ships not utilizing OTSR, when making a transoceanic voyage in the 
Northern Hemisphere, usually travel a more northerly route in the summer 
months and then shift to a southerly route in the stormy winter months. 
This usually results in avoiding the largest storms but at a great cost 
in time and money. Reduction in weather damage is not the whole solu- 
tion, however, since the ideal situation is to save time and increase 
the distance traveled per day. Often times, the weather situation may 
be such that the off-season route has the best sea conditions. 

At the present time a 5-day forecast of the surface weather phenom- 
enon is the maximum utilized in ship routing. If this could be extended 
say to 15-20 days, an entire voyage could be routed across the Pacific 
Ocean to avoid storm areas and reduce transient time. Operationally, 
the procedures for ship routing are slightly different in the Atlantic 
and Pacifie operating areas. More emphasis is placed on storm track 
trends and climatic means in the Pacific since the 5-day retdtart covers 
less than half the normal crossing. 

It is the purpose of this paper to derive extended range forecasts 
of wave heights, directions and periods for 15-20 days in advance of a 
proposed sailing across the Pacific Ocean. 

Naturally, constant surveillance would have to be maintained over 
the routed vessels as the 5-day forecasts are up-dated and when sea 
conditions warrant, the routes would be modified to indicate a superior 


route. 


2. COLLECTION OF DATA 

The U.S. Weather Bureau publishes a 5-day forecast of the mean sea 
level pressures around the Northern Hemisphere and this is presented on 
individual maps for each of the five days. These maps are distributed 
on Monday, Wednesday, and Friday of each week and are portrayed on a 
polar stereographic 1:50,000,000 scale. The forecasts are based on a 
mean circulation method by averaging ten twice-daily values of pressure 
at standard intersections of latitude and longitude on synoptic sea 
level weather maps. (1 | the Weather Bureau also makes 30-day forecasts 
using methods similar to those employed in preparing the 5-day fore- 
casts. 

These maps were utilized to obtain a surface wave forecast for an 
extended period of time across the Pacific Ocean. The period of time 
selected was 26 July 1966 to 06 August 1966, a total of 12 days. The 
U.S. Navy Fleet Numerical Weather Facility located at Monterey, Califor- 
nia publishes 24 and 48-hour wind-wave prognoses which were used for the 
first two days of the forecast, 26 and 27 July 1966. [2] Bor 28.29 , .agd 
30 July 1966 the last three maps of the U.S. Weather Bureau 5-day fore- 
caSt series were applied. For the 3lst of July through 6 August the 
30-day mean sea level prognostic chart dated 15 July 1966 was utilized 
on a day to day basis. 

In the following explanation, focus is placed only on one map in 
showing how the surface wave forecast was derived. The remaining days 
of the forecast were treated in a Similar manner. 

The first step in arriving at the forecast was to use the FNWF 
FAREB program. This program was written by Harry N. Farnsworth on 1 


May 1965 at the Fleet Numerical Weather Facility. The basic purpose 


of this program is to reconstruct any scalar field in memory in a CDC 
1604 computer as a 3969-point field. An attempt will not be made to 
explain each step of this program in detail, however, input and output 
values along with an example for preparation of a field will be dis- 
cussed. 

Input values were extracted from the source chart by employing a 
plastic overlay of the standard FNWF 63 x 63 grid. If a grid overlay 
were not available latitude-longitude input could be used. A proper 
interval should be selected in order to combine sufficient pattern de- 
lineation with a minimum of data inputs. The initial guess field is 
considered in the interval selection and since a guess field of zero 
was used, an interval of two grid points or every other point was chosen. 
Then by superimposing the grid overlay on the source chart an initial 
grid point was established and grid point values were read along a J-row 
at a specified point spacing (every other point). A new initial point 
may be established whenever desired to reorient the input data. The 
input points selected were denoted by dots and dashes on the FNWF over- 
lay. Since the source chart was in millibars, the values were read to 
the nearest millibar interpolating where necessary. When all values had 
been read from the source chart in the desired area and tabulated they 
were converted to 1000-mb heights using a standard conversion of 7.5 mb 
per 200 feet with a reference point of 1000 mb equal to zero feet. All 
values had to be coded as four digit numbers for entry on the standard 
IBM card. The values were assumed to be positive unless preceded by a 
minus sign. There were some low centers present with central pressures 
below 1000 mb and these were coded accordingly. 


In order to assemble the card deck the identification and scaling 


information was prepared. On the initial card the titling information 
and initial guess value were placed. The standard procedure is for the 
titling information to consist of the date-time group, some identifying 
number, in this case D1000, and then any explanatory information that 
might be desired such as name, etc. 

The program was then put into operation and the final field con- 
structed through smoothing, relaxation, and a combination of the tw. 
A grid print in feet of the height of the 1000-mb surface was produced 
and the field packed and written on magnetic tape. As mentioned pre- 
viously the main objective of this program was to get the data on tape 


in an organized field in order to go on to the next Step. 


3. WIND EXTRACTION 

The next step in arriving at the wave forecast was to utilize the 
FNWF wind extraction program. By using the 1000-mb constant pressure 
surface as the input this program computes the geostrophic u, v wind com- 
ponents, as well as the direction and speed of the wind by solution of 
the geostrophic wind equation. There are, however, some limitations to 
this program. Wind speeds that are computed by the geostrophic wind 
equation are generally greater than actual speeds in regions of cyclonic 
trajectory curvature and less than actual speeds in regions of anticy- 
clonic curvature. 

This program was run on the 3200 computer at the Fleet Numerical 
Weather Facility. An example of a 3200 run sheet is shown in Appendix 
I. Some modifications had to be made in order to read the mean D1000 
field and the climatological D1000 field taken from the U.S. Weather 
Bureau records. These modifications were made on a 3200 octal patch 
and an example of this is shown in Appendix II. The output of this 
program was packed and written on magnetic tape for an input to the wind 


waves program. 


10 


4, WIND WAVES 

The next step was to derive the final product. The Fleet Numerical 
Weather Facility Wind Waves program with slight modifications was run on 
the Control Data Corporation 1604 computer to accomplish this purpose. 
In order to eliminate some of the steps not needed in this program a 
flex-tape patch (Appendix III) was made and inserted into the computer 
at the appropriate time. 

The program input consisted of three magnetic tapes. Tape number 
one contained the 1000-mb field obtained in the FAREB program previously 
mentioned. Tape number two contained the climatological sea surface 
temperatures for the month of July. Tape number three was comprised of 
the u and v components of the wind. A fourth tape, blank pool tape, 
was used in order to copy the output of the program. This was necessary 
in making a difference field with the actual analyzed waves in order to 
determine the merits of the 5 and 30-day forecasts for preparing long 
range wave prognoses. 

The Waves program corrects the geostrophic wind directions for 
surface friction by assuming a constant cross-contour flow of 15 de- 
grees toward lower pressure. A stability correction is then made by 
reducing the wind speed depending upon the general air mass type and 
the sign of the surface temperature advection. For example, a mass of 
Maritime tropical air blowing toward colder sea surface temperatures 
would warrant a greater reduction in wind speed than a mass of maritime 
polar air blowing toward warmer sea surface temperatures. The duration 
of the wind is determined at each grid point by computing the length of 
time the wind has been blowing from the same direction (within * 22 de- 


grees). The maximum duration has been set at 18 hours according to the 


ll 


French school of thought which considers all seas to be fully developed 
after this length of time. An effective wind speed is then obtained by 
an averaging technique. The significant wave height and wave period is 
computed using two singular relationships established at the Fleet Nu- 


merical Weather Facility: 


= x. v2 
Hy 73 = K,V4D + KV (4.1) 


11/3 V(K3 52 K),D) ate Ke (4.2) 


where Hy /3 and T1/3 are the significant wave height and period, respec- 
tively; V is the effective geostrophic wind speed; D is the duration 
time and K, to Ks are empirical constants tuned to reduce the root mean 
square error. 

An example of a 1604 run sheet is shown in Appendix IV. The tape 


set-up and instructions listed were actually used to run this program. 
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5. VERIFICATION 

It was decided at the onset that once the forecast wave pattern was 
derived a comparison would be made of the derived field and the actual 
analyzed waves for the specified period of time. To accomplish this 
purpose an A-B program, called SAK and written by R. McDonald at Fleet 
Numerical in 1962, was utilized. This program computes a difference 
field using any two given fields, A and B, a pillow and root mean square 


error for the entire grid. Pillow and RMSE are defined as follows: 


Pillow 


x 
ie ee x = Dependent on number (5.1) 
nat of points differenced 


x 


(52) 


x 
2 
RMSE = E [cam = Pillow | 
n 
n=l 


x 





A further comparison was made of the actual analyzed waves and a 
wave forecast based solely on climatology. Sea level pressure normals 
were obtained from the U.S. Weather Bureau and a D1000 field was de- 
rived. From this field a wave forecast was made following the same 


basic steps as explained above. 
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6. CONCLUSIONS AND ACKNOWLEDGEMENTS 

Appendix V portrays the results of the verification that was out- 
lined in Section 5. It is readily apparent that the waves derived from 
the 5 and 30-day forecasts show more skill than those waves derived 
solely from climatological data. The reason the values of RMSE in both 
cases might seem to be large is that data was extracted only in the 
Pacific Ocean region and the rest of the grid left blank. When the SAK 
program computes this error the whole grid is taken into account result- 
ing in the large values. Based on this fact, it is felt that the pro- 
cedures outlined in this paper can be used in extending wave forecasts 
beyond those that exist today with a good degree of certainty. 

The author wishes to thank Professor George J. Haltiner, Chairman, 
Department of Meteorology and Oceanography, U.S. Naval Postgraduate 
School, Monterey, California for his guidance and a special thanks goes 
to Mr. Norman Stevenson, Data Dynamics, Monterey, California for his 
invaluable assistance in this endeavor. 

The author also desires to acknowledge the cooperation of Dr. 
Jerome Namias and Mr. James O'Connor of the Extended Range Forecasting 
Section of the U.S. Weather Bureau in providing the long-range weather 
forecasts used in this study. 

Finally, the author wishes to thank the personnel at Fleet Numeri- 
cal Weather Facility, Monterey, California for the use of their compu- 


ters and for providing the data. 
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APPENDIX I 


3200 Run Sheet 


Run # 12 Date 9-26-66 Est. Time 10 Stop if exceeds 
Programmer Truax Actual Time 





Rush Rout ine Rockpile Listing Only 


C ) BPI IN OUT List & Binary Deck 
or) BPI IN OUT List & Mini Tape 
( ) BPI IN OUT Stack for Transmit 


() I Fs a Program Requires: Drum Printer 


1. Load: Binary Tape Binary Deck Patch Cards 


2. Load; Drum Records 
Tape Lib Records 


3. I&C Loc aoe = 


Loc = ; Loc = 





4. Set JK : SK_: Other 


5. Set Wds 13,14(Dtg) Wd15_ Wdl6  Wdl7__ 


Set A 


ul 
O 
iT 


Bl = B2 B3 = 


BPI = BPO = STO = 


1) 
ct 
wy 
4 
cr 
a) 
uN 


Expected Final Dump Yes No 


Other Instructions: 


(1) Set up for "FWBXB" using tape H-24 (200 Density) as input. 
(2) Read in Bi-Octal correction cards. 

(3) Run wind 31A on output fields from FWBXB. 

(4) Save output on enclosed ST-500 tape. 
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02435 


02436 


02437 


02440 


02441 


02442 


02443 


O2444 


02445 


02446 


02447 


02450 


02451 


02452 


02453 


02454 


01007 


APPENDIX IIL 


64446561 


45202020 


00000170 


000000 36 


00567000 


00160400 


63121220 


11005305 


646 34371 


44612320 


00000170 


000000 36 


00567000 


00160400 


63121220 


11005 305 


00010201 


3200 OCTAL PATCH FOR WIND EXTRACTION 
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00705 


00706 


00707 


00710 


O07 11 


00712 


00713 


00714 


G07 1a 


00632 


00751 


01056 


01057 


01060 


01061 


01062 


01063 


01173 


01175 


ODL 


OZ 


O1S20 


01724 


02007 


ooooo;qcoooorcorc;cveoco cr corr oorrroaoooooce°orooorFrooCcoeorfroooFrFC]°o 


1604 FLEX PATCH FOR WIND WAVES 


12,35:32 
Japaayeyh 
05150 
01400 
22351 
32170 
05150 
01400 
62170 
42007 
05150 
00100 
61445 
71264 
05150 
00200 
00000 
00000 
00746 
00000 
01053 
00000 
591626 
T2552 
07 6 34 
01056 
00004 
00000 
S2170 
00000 
07600 
01061 
01157 
00000 
O22 7 
12582 
C1276 
2532 
01315 
00000 
P2552 
22366 
01720 
00000 
02564 
00000 
02564 
00000 


APPENDIX III 


02602 


02566 


02600 


03553 


01026 


01027 


010 30 


01031 


01032 


01033 


01034 


01035 


01036 


01037 


01040 


01041 


01042 


01043 


01044 


01045 


01046 


00716 


00717 


00720 
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ooo o°0c co 0coc co 0c C0 0C O00 0O0 O0 COFOWOnNnNOOUVVOYFP FwwohOreOoFoeo0o0C OF COFCO O 


03501 
00000 
00000 
12332 
01763 
[2332 
01211 
00000 
10041 
T2002 
05150 
01400 
00000 
00000 
00000 
00000 
10041 
00000 
00000 
00014 
00006 
T2532 
00001 
00000 
00003 
01033 
Oule/o3/ 
Oe 32 
01041 
00000 
22328 
010 34 
01044 
01040 
01026 
00000 
2 332 
01034 
01041 
01040 
00716 
00000 
00745 
007 37 
00745 
007 33 
2170 
007 35 


00721 


00722 


00723 


00724 


00725 


00726 


00727 


007 30 


007 31 


00745 


01214 


01215 


01216 


OP217 


03560 


00761 
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APPENDIX III (CONT'D) 


10041 
12514 
05150 
00100 
00733 
00000 
61445 
007 35 
10041 
12514 
05150 
00200 
00733 
00000 
32170 
00735 
00627 
00000 
DO7 33 
00000 
00000 
00000 
00000 
00000 
00000 
00000 
01026 
00000 
00001 
03554 
00000 
01270 
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APPENDIX IV 


1604 Run Sheet 


Run #/ Date 9-16-66 Est. Time 10 Stop if exceeds 


Programmer Truax Actual Time 
Rush Routine X  Rockpile Assembly: AR Scrap 
Ch 3/4 Bioct. Bin.Cards 
Ring } 


(1) 1000 mb BPI 200 in. out X 


(2) SST BPI 200 in. out xX en 5/6 


Ring 


( ) BPI ine out (3) Pool BPI 556 in. X out 


BPI in. out (4) Ujmeeset 200 in our: 


1. Load Drum library records "WAVES" (Operational program) 


2. Load bioctal tape(s) Flex tape(s) one 


3. I&C Location ACL) ACR) 





A(L) A(R) 
4. Set Wd. 10 = Wd. 16 Pos __Neg 
5. Set Breakpoint _ JK_1 2 


6. Set A — ae __ sect 0 





7. Start P = 1026 Expected Final 616 Dump Yes X_ No 


8. Dump 600 to 10000 Instr/Data 


Other Instructions: 


1. Start at 1026 program will stop at 616, hit go 
2. Save pool tape if program runs 
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Date 


Zo J why 
Zo July 
sO July 
peo uly 
1 August 
2 August 
3 August 
4 August 
5 August 


6 August 


APPENDIX V 


COMPARISON OF RMSE & PILLOW 


Derived Waves 


RMSE (FT) 
2 


3 


PILLOW (FT) 


0 


il 


Ze 


Climatological Waves 


RMSE (FT) 
6 
6 


6 


PILLOW (CFT) 
3 
3 


3 
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